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Primordial black hole formation and hybrid inflation 

David H. Lyth 
Department of Physics, Lancaster University, Lancaster LAI 4YB, UK 

We calculate the contribution to the curvature perturbation, that is generated while the waterfall 
field grows exponentially (which can occur only if the waterfall field mass ism> H). We consider 
the upper bound on its spectrum coming from constraints on the abundance of primordial black 
holes, taking account of possible non-gaussianity. The constraint is satisfied, but extrapolation of 
our result to the regime m ~ H suggests that the constraint might not be satisfied there. 



I. INTRODUCTION 



On cosmological scales, the primordial curvature perturbation Q is very small with spec- 
trum 'Pf(fc) — (5 x 10~ 5 ) 2 . But V^k) may rise as the wavenumber k increases, to some 
much bigger peak value 7\(fc pC ak)- If that happens, primordial black holes may form as 
fcpoak enters the horizon, and cosmological bounds on their abundance translate to an upper 
bound on ^(fcpoak)- 

The bound is usually discussed under the assumption that C is nearly gaussian, giving 
\mJ • very roughly ^(fcpcak) ^S 10 -2 . In the first part of this paper we discuss the bound assuming 

r~| | instead that on scales around fc pea k 



C(x)=±( 3 2 (x)-(. 9 2 )), (1) 



9*' 
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with g gaussian. We show that the bound becomes very roughly Vq (fcpeak) ^ 1CT 3 for the 
positive sign and 7\(fc P cak) ^S 1 for the negative sign. 

It is known |l| that the form ([1]) can hold if the rise in V^k) is generated after inflation 
by a curvaton-type mechanism. For the curvaton mechanism itself the sign is positive but 
it could be negative more generally. With a curvaton-type mechanism, fc pca k corresponds to 
the horizon scale when the mechanism ceases to operate and black hole formation follows 
00 . immediately. 

A different possibility is for the rise to be generated by the waterfall held of hybrid 
inflation. A calculation has been done for a particular case 2], where simple expressions 
were found for £ and Vq. In that case £ is gaussian and the black hole bound is well satisfied. 
The main purpose of this paper is to generalise the calculation of [2] , to cover essentially 
any hybrid inflation model permitting exponential growth of the waterfall field. 



II. BLACK HOLE BOUND 



5_j 

The bound that we are going to consider rests on the validity of the following statement: 

if, at any epoch after inflation, there are roughly spherical and horizon-sized regions with 
£ significantly bigger than 1, a significant fraction of them will collapse to form roughly 
horizon-sized black holes.* 1 The validity is suggested by the following argument: the over- 
density at horizon entry is 5p/p ~ £, and if it is of order 1 then 5p ~ p = 3MpH 2 . The 
excess energy within the Hubble distance H~ x is then M ~ H~ 3 p ~ Mp/H, which means 
that the Hubble distance corresponds roughly to the Schwarzchild radius of a black hole with 
mass M . The validity is confirmed by detailed calculation using several different approaces, 
as summarised for instance in [3|. 

Before continuing we mention the following caveat. Practically all of the literature, as 
well as the simple argument just given, assumes that £ within the region is not very much 
bigger than 1. Then the spatial geometry within the region is not too strongly distorted and 
the size of the black hole is indeed roughly that of the horizon. In the opposite case, the 



* 1 As in Eq. |[TJ we are choosing the background scale factor a(t) so that the perturbation C, = <5(ln a(x,t)) 
has zero spatial average. 



background geometry is strongly distorted and the wavenumber k defined in the background 
no longer specifies the physical size of the region at the epoch aH = k of horizon entry [4| . 
An entirely different discussion would then be necessary, which has not been given in the 
literature. As the opposite case does not arise in typical early-universe scenarios we ignore 
it. 

We are interested in the case that V^(k) has a peak at some value fc pea k, and we assume 
that the width of the peak in lnfc is roughly of order f so that (£ 2 ) ~ ^(fcpoak)- We will 
focus on the case that £ is generated by the waterfall field perturbation. Then V^(k) ex fc 3 
below the peak with an exponential fall-off above the peak, giving^ 2 

(C 2 ) ~ / V c (k)dk/k = -2>c(W), (2) 

but the numerical factor is not very important. What matters is that regions with £ > 1 
that might form black holes will be rare if Pf (fc pea k) is sufficiently far below f , 

Observation demands that the regions must indeed be rare, because it places a strong 
upper bound on the fraction of f3 of space that can collapse to form black holes, on the 
assumption that the collapse takes place at a single epoch as is the case in our scenario. A 
recent investigation of the bound is given in [3j], with extensive references to the literature. 
The bound can arise from many different physical effects, depending on the epoch of collapse, 
and is subject to many uncertainties. Typical bounds are roughly in the range 

1CT 20 < /3 max < KT 5 , (3) 

and we shall take that to be the range in what follows. To bound Pf (fcp ea k), we shall require 
/ < /?max, where / is the fraction of space with ( > 1. 

That fraction can be calculated from (£ 2 ) if we know the probability distribution of C( x )- 
The standard assumption is that it is gaussian. Then 

/ = crfc(f/V2(C 2 ) 1/2 ), (4) 

and using the large-x approximation erfc(x) ~ e~ x I 2 we find (£ 2 ) ~ l/21n(l//). For the 
range © this gives (with Eq. ©) T c (k pca . k ) < O.Of to 0.04. 

If instead C has the non-gaussian form ([I]) with the plus sign we have 



(C 2 )-2(.g 2 ) 2 



1 



(5) 



.21n(l//) 

which with Eq. © gives 7\(fc peak ) <6x I0~ 4 to 5 x I0~ 3 . 

The situation when £ has the form ([I} with the minus sign is quite different. There is 
now no region of space where £ > (g 2 ), and / <C 1 now implies some bound (g 2 ) -1<1 
which is practically equivalent to (g 2 ) < 1. With Eq. ([2]) this corresponds to Vq (fc pe ak) < 6. 

III. EVOLUTION OF THE WATERFALL FIELD 

We adopt the notation and basic approach of [2]. In the first two subsections these are 
summarised, and then we present our new calculation of evolution of the waterfall field %. 

A. Standard hybrid inflation 

We consider the usual hybrid inflation potential, 

Vfax) - V + V(<t>) + ±m 2 (cl>)x 2 + \\x 4 (6) 

m 2 (<P(t)) = g 2 4> 2 (t) -m 2 ^ g 2 (cb 2 (t) - 2 ) , (7) 



# 2 If the integral fails to converge at small k we impose a cutoff corresponding to the size of the observable 
universe. The contribution from small k then has a negligible effect. 



with < A <C 1 and < g <C 1. We consider the standard case m ^> H. The effective 
mass-squared m 2 ((j)(t)) of the waterfall field goes negative when the inflaton <fi falls below 
<t> c = mj g. Then \ moves towards its vev, marking the beginning of what is called the 
waterfall. The waterfall is deemed to end when \ approaches its vev, which usually marks 
the end of inflation. 

At least after the observable universe leaves the horizon, the inflaton potential V(4>(i)) 
is supposed to have positive slope, and the inflaton is supposed to have zero vev so that 
V(0) = 0. Also, it is assumed that V(<p) < V . 

The requirements that V and dV/dx vanish in the vacuum give the vev xo and the 
inflation scale V ~ 3M^H 2 : 

x g = ^~12JW|fl 3 /m a , V = ^c3M*H 2 . (8) 

The waterfall is supposed to begin with an era during which the evolution of x is linear. 
Choosing a gauge whose slicing corresponds to uniform </>, 

Xk + 3#Xk + [{k/af + m 2 (<p(t))] Xk = 0. (9) 

The energy density and pressure of x are 

Px = ™\<P)X 2 + \X 2 + \\V X \ 2 (10) 

p x = - m \4>)x 2 + \x 2 + \\V X \ 2 - (11) 

In [2j we considered a regime of parameter space is which (i) the waterfall takes much less 
than a Hubble time, (ii) m 2 (t) oc t during the linear era and (hi) x is growing exponentially 
by the end of the linear era. We found simple formulas for the contribution £ x of x to the 
curvature perturbation to at the end of the linear era, and for its spectrum V^ ■ In this 
paper show that very similar formulas hold if we assume only exponential growth. 



B. Exponential growth 

Since inflation by definition continues during the the waterfall, H will not vary much 
and we set it equal to a constant to simplify the presentation. Then conformal time is 
r] = —1/aH and |J9]) can be written 



! ( a Xk) . ,j 



Wfc«Xk = 0, (12) 



dr] 2 
with 

uj 2 k {r}) = k 2 + a 2 rh 2 (t), m 2 = m 2 (t) - 2H 2 , m 2 (t) = g 2 cf> 2 (t) - m 2 . (13) 

For sufflcientlly small k, we can set (Jj. ~ oj 2 =0 — a 2 rh 2 . Then w^ switches from positive 
to negative before <f> = <p c , but presumably not long before since m 3> H . For k 2 > the 
switch is later. For the scales that we need to consider, we assume that there are eras both 
before and after the switch when tu 2 , satisfies the adiabaticity condition d|w 2 |/c?7/ <C |w^|. 

During the adiabaticity era before the switch we take the mode function to be 

a Xk c- (2Lo k (r]))- 1/2 exp (-i / uj^dr]) , (14) 

which defines the vacuum state. During the adiabaticy era after the switch 

a Xk - (2K(r,)|)- 1 /2 exp ( /"' 1^(77)1^ ] , (15) 



where the subscript 1 denotes the beginning of the adiabatic era. The displayed prefactor 
holds [2j only if m 2 (t) ex t and H(t — ti) <C 1 but its precise form doesn't matter. All we 
need is for it to vary sufficiently slowly that the growth is dominated by the exponential. 
We call the era during which that is true the growth era. 

During the growth era, the adiabaticity condition is equivalent to the three conditions 



m(t) ~ Tn(t), 
2H 



\m(t)\ 

d\m{t)\ 



\m{t)\ 2 dt 



< 



< 



1 



o(t)|m(*)| 



o(t)|m(*)| 



3/2 



3/2 



(16) 
(17) 

(18) 



The adiabatic era begins when all three conditions are first satisfied. 

By virtue of Eqs. ([TBI and (fl"7|). we have from Eq. ([TBI Xk — \m{t)\xk- Except near places 
where x(x, t) vanishes, this implies 



X(x,t)~|m(t)| x (x,t).. 
In the regime k <C a|m(t)| we have 

1 k 2 



\u k \ ~ a\m(t)\ 1 



2a 2 \m{t)\< 



(19) 
(20) 



giving 



Xk (t) ~ X k=o(t)e- k ^ 2k ^\ Xk=0 (t) ~ (2a 3 |m(t)r 1 / 2 exp (J dt\m(t)\) , (21) 



where 



kl[t) = 



dt 



tl a 2 \m(t)\ 



(22) 



During the growth era, Xk(0 is classical and proportional to the mode function Xk(t)- 
The same holds for x(x, i), except near places where it vanishes. The spectrum is P x — x\- 
Also, 



(x 2 (t) 



47T 



r-P x (o,r) 



dkk 2 



-(k 2 /kl(t)) 



(27r)- 3 / 2 P x (0,r)^(t), 



and using the convolution theorem we have [2] for k -C fc, 



#:s 



1 ' ,2/ + N\2 r , /, , m3 



n x2 (t,fc) = — ( X 2 (i)) [fe/fc*(r)] d , 



(23) 



(24) 



with Vs x ^{t, k) falling exponentially at fc 3> fe*. 

By virtue of Eq. (fT8"|). the change in |m(t)| in time |m(t)| _1 is negligible and so is the 
change in a. Setting t = t\ + |m(i)| _1 we get k 2 (t) ~ a 2 (i)|m,(i)| 2 . But at this epoch the 
growth era has hardly begun, and subsequently k*(t) decreases while a\m(t)\ increases. So 
we really have k 2 (t) -C a 2 \m(t)\ 2 during the growth era, and Xfe falls exponentially in the 
regime k*{t) < k < a\m(t)\ which means that fc*(t) is the dominant mode. 

We are mostly interested in an epoch just before the end of the linear era, which we 
denote by a subscript n i- We denote fc*(i n l) simply by k*. Let us define N n \ = H(t n \ — t\). 
If AT nl < 1, 



1< 



\m(h)\ 
H 



< 



a(t nl )H 



< 



\m(t n \)\ < m 
H ~ H 



(N nl < 1). 



(25) 



* 3 As usual P x = (2tt 2 J k 3 )V x , with both P and V referred to as the spectrum. 



If instead N n \ > 1, fc*(£) levels out after H(t — t{) ~ 1. We therefore have, whatever the 
value of iV n i, 

Since |m(i)| 3> H, the upper bound on fc* imples (fc*/a)xfc <C Xfe = l m (0IXfc- The spatial 
gradient of x( x j wu l therefore be small compared with the time-derivative, except near 
places where x = 0- Hence 

Px ~ -l| TO (t)|V + ^ 2 ^0 (27) 

p x ~ i|m(t)| 2 X 2 + ix 2 ^ |m(t)|V * X 2 . (28) 

To evaluate p x we can use the local continuity equation (valid because the spatial gradient 
is negligible), 

p x = -3H(p x +p x )~-3Hp x , (29) 

Px ~ -(3#/2|m(t)|)p x . (30) 



giving 



IV. CURVATURE PERTURBATION 

Now we calculate the contribution of x to the curvature perturbation, following closely the 
procedure of [2J]. The curvature perturbation is given by ((x,t) = 5N(x,t), where N is the 
number of e-folds of expansion from any initial slice with a(x, t) = a(t) (flat slice) to a slice 
of uniform p at time t. The contribution of x during the linear era is Cx( x ) = $N(yi, t n \, ii), 
where N(x,t,ti) is the expansion from a slice of uniform p just after the beginning of the 
linear era, to a slice of uniform p just before its end. We are working in a gauge where 
5(f) = so that p(x, t) = p x (x,t) + p$(t). Since \Sp x (x,ti)\ <C \5p x (x,t n i)\ we have 



r 5p x (x,t n i) _ 1 8p x (x,t n \) 

Kt) Z{x 2 {tm))+HtY 



Cx(X; t) = _ g "*Wni; = x _^x^u_ ; (31) 



where p(i) is the spatial average of p(x, £). 

Using the equations of the previous section this gives for k <C fc* 

> r „s g (x 2 fa0) *x 2 (x,tm) . . 

^ J 2|m(t n i)|& a (tni)> + 2 (tm) (X 2 (im)) ' l j 



and 



*>Cx(*) 



H (X 2 (im)) 



2|m(t nl )|( X 2( inl )) + 02 (tnl) 



fc 



3 



(33) 



At fc 3> fc*, 7^ is negligible because Xfc is. The spectrum therefore peaks at fc ~ fc*. 

If (x 2 (^ni)) 3> <^ (tnl) we have p ~ p x , which means that the slice of uniform p is practically 
the same as the one of uniform p x . In turn, that is the same as the slice of uniform %. 

In general the formula £ = SN holds only after smoothing p and p on a scale big enough 
that the local continuity equation is satisfied. In our case though, that equation is satisfied 
on all of the scales fc 3> fc* on which £ x is significant. Therefore, the formula Cx( x ) = 
6N(x, t n i,ti) makes sense on all of these scales. In the opposite regime fc <C fc*, ( x will be 
negligible simply because <5x 2 is, its exponential growth not having begun. Thefore, V( 
peaks at fc* whether that scale is super-horizon or sub-horizon. 

We want to see whether the black hole bound of Section [TT] is satisfied by Eq. (|3"3"|) . The 
situation is simple if fc* is super-horizon. Then, £ is of the form Eq. {T]) with the minus sign, 
and the black hole bound is Vc, (fc*) < 6 which is well satisfied. 



If instead k* is sub-horizon, we have to remember that the black hole bound refers to 
horizon-sized regions. To apply it, we must drop sub-horizon modes of (, x . Estimating the 
bispectrum, trispectrum as in |2|, one sees that this makes C, x nearly gaussian. The black 
hole bound is therefore roughly V$ X (H) < 10 -2 . Since fc* >• fc C nd and \m(t)\ ^> H, it too is 
presumably satisfied. 

Earlier calculations of £ x are reviewed in [2[ , including the one in Q . A continuation of Q 
has since appeared |6j. These papers consider the potential V(<j)) oc 4> 2 , for which m 2 (t) oc t 
is a good approximation. Their expression for ( x (Eq. A(10) of @ and Eq. (6.15) of |6|) is 
the same as our Eq. (|32l) . except that the middle factor is missing. However, the meaning 
of Sx and (x 2 ) in the final factor is different in their expression from that in Eq. ([512]) ; \ m 
their case is smoothed on the horizon scale whereas in our case it is not. 



V. CONCLUSION 

We have calculated the contribution to £ generated during the linear era of the waterfall, 
assuming that such an era exists and that exponential growth takes place during that era. 
Such growth can occur only if the tachyonic mass of the waterfall field is m 3> H (standard 
hybrid inflation) which we therefore demand. 

The calculation generalises [2|, to the case of an arbitrary time-dependence of m 2 (t) 
subject to it being slow enough to allow exponential growth. Because the dependence is 
arbitrary, we cannot repeat the detailed investigation given in |2j, of the parameter space 
within which our calculation will apply. That could only be done with a specific inflaton 
potential V (</)), which among other things would determine m 2 {t). 

We have, as usual, considered only super-horizon modes of the curvature perturbation 
£. The usual super-horizon treatment of £, based on the local energy continuity equation 
does still hold for sub-horizon modes during the waterfall. But we have not used these 
modes, and to do so would require evolution of the cosmological perturbations through to 
the post-inflation era. In particular, that would be necessary if we were to discuss the 
possible formation of black holes whose size is much smaller than the horizon scale at the 
end of inflation. 
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